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Abstract. Correction due to finite speed of light is among the most inconsistent 
ones in absolute gravimetry. Formulas reported by different authors yield corrections 
scattered up to 8 /iGal with no obvious reasons. The problem, though noted before, has 
never been studied, and nowadays the correction is rather postulated than rigorously 
proven. In this paper we make an attempt to revise the subject. Like other authors, 
we use physical models based on signal delays and the Doppler effect, however, in 
implementing the models we additionally introduce two scales of time associated with 
moving and resting reflectors, derive a set of rules to switch between the scales, and 
establish the equivalence of trajectory distortions as obtained from cither time delay or 
distance progression. The obtained results enabled us to produce accurate correction 
formulas for different types of instruments, and to explain the differences in the results 
obtained by other authors. We found that the correction derived from the Doppler 
effect is accountable only for | of the total correction due to finite speed of light, if no 
signal delays are considered. Another major source of inconsistency was found in the 
tacit use of simplified trajectory models. 
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1. Introduction 

Modern absolute gravimeters measure gravity acceleration by laser tracking and 
subsequent analysis of the trajectory of the free falling test body. The laser beam 
(L) (fig. 1) split by the mirror (BS) and reflected simultaneously from the fixed (FCC) 
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Figure 1. Functional diagram of an absolute gravimeter 

and the free-falling (FFCC) corner cube qua test body, creates the interference signal 
converted by the photo detector (PD) into the fringe wave (FW), processed by the 
unit (PU) to obtain a set of time-distance coordinates of the falling body and calculate 
the gravity acceleration. The set of coordinates {S^Tj}, also called "levels", makes up 
the measurement schema. If the falling body is tracked on both upward and downward 
parts of the trajectory, the schema is rise-and-fall one, otherwise the schema is free-fall. 
There are schemas with two, three, or four levels, as well as multi-level schemas, in 
which the levels may be equally spaced in time (EST), or in distance (ESD). Although 
the majority of absolute gravimeters now use multi-level schemas, in this paper we 
consider other schema types as well. One reason for this is that older schemas are still 
used, albeit rarely [1]. Another reason is that theoretical results for non-multi-level 
schemas are often extended to the multi-level case [2, 3, 4], not realizing the fact and 
its consequences. 

Due to finite speed of light, the laser beam reflected from the free-falling corner 
cube, delivers the information of its position with some delay, resulting in a distortion 
of the measured acceleration, which therefore requires a correction. Though the models 
of this phenomenon are pretty straightforward, the correction due to finite speed of light 
is perhaps the most controversial one in the theory of absolute gravimeters. There are 
more results published on this correction, than on any one else, that have later been 
disputed or amended by other authors. Generally, corrections are known to depend on 
the measurement schema of the gravimeter. For the finite speed of light, however, the 
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Table 1. The values of k in formula (1) obtained by different authors. 

k author 

| Murata [2] 

1 Tsubokawa T. [4] 

| Hanada [4] 

| Kuroda & Mio [5] 



obtained corrections differ significantly even for the same schema types. Let's examine 
some results published for EST schemas. Neglecting the correction sign and the test 
body's initial velocity! , all authors agree that the correction is given by the formula 

Ag c = k 9 -^, (1) 
c 

where g is gravity acceleration, c is speed of light, T is the duration of the trajectory, 
k is the disagreement factor. The values of k obtained by different authors are shown 
in the table 1. As the duration of the trajectory T in modern gravimeters is 0.1. ..0.2 s, 
the values of k from the table yield the corrections different as much as 5.5 /zGal § ||. 
This uncertainty exceeds significantly the precision level of modern instruments, raising 
the question of the "correct" correction. 

Another way of accounting for the finite speed of light is to modify the measured 
time intervals Tj. The gravity acceleration is obtained by the model [6]: 

S i = z + V Q f i + g Q f?/2, (2) 

where 

f i = T i -^^. (3) 

c 

Here z , V , and g stand for initial coordinate, velocity, and acceleration of the test body 
respectively. This model is non-linear by the parameter zq, independent determination 
of which is rather problematic. It's also unclear if the gravity acceleration obtained with 
the above model can be matched to any value from the table 1. 

In this paper we consider several approaches to the correction due to finite speed 
of light and suggest consistent ways for its implementation. In doing so, we uncover 
the reasons of the divergence of the results obtained by other authors. The paper 
has following structure. In chapter 2 we discuss the approach we use to account for 
disturbances in the corner cube trajectory. In chapter 3 we consider several models 
of influence of the finite speed of light on the measured gravity value and establish 
connections between the models. In chapter 4 we find and investigate the corrections 

| The sign of the correction depends on the position of the interferometer with respect to the test 
body. The initial velocity was considered not by all authors, so we neglect it for this example. 
§ 1 ,uGal = 1(T 8 m/s 2 

|| The result of Hanada [4] includes a miscalculation, which, if corrected, makes his k equal to |, 
increasing this difference to 7.5 /xGal. The result of Hanada is discussed in chapter 5.6. 
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for different measurement schemas. The paper is concluded by the chapter 5 where 
we review the results obtained by other authors. Elements of the theory of absolute 
gravimeters that require additional comments for this paper, have been carried into the 
Appendices. 



2. Accounting for disturbances in the test body motion 

We consider the measured gravity value as weighted average of the test body acceleration 
9it) [7]: 



l 

Jg(t)w(t)dt, (4) 



i 

I 



T 

9 

o 

where g is the measured gravity, g(t) is the acceleration of the test body, T is the duration 
of the trajectory, w(t) is the gravimeter's weighting function, which is an analog of the 
impulse response function of a linear system. The intrinsic property of the gravimeter's 
weighting function is its unit square: 

T 

w(t) dt = 1. (5) 
o 

If along with its normal value go, the acceleration of the test body contains a disturbance 
Ag(t), the measured gravity, according to (4), will be 

T T 

g = J \g + Ag{t)) w(t) dt = g + J Ag(t) w(t) dt. (6) 
o o 
So, if the gravimeter's weighting function is known, the additional acceleration measured 
due to the disturbances Ag(t) can be found as 

T 

Ag~{t) = J Ag(t)w(t)dt. (7) 
o 

Taken with the opposite sign, the above expression yields the correction for the 
considered disturbance. 

If the disturbance is expanded into the power series like 

Ag(t) = J2^nt n , (8) 

the accounting for it is greatly simplified. In this case, the measured additional 
acceleration, according to (7), is 

T 

Agjf)= J (5>n* n )™(0d* = 5>nC n , (9) 
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where 

T 

C n = Jt n w(t)dt, n = 0,l,2,... (10) 
o 

is the ra's averaging coefficient of the gravimeter. These coefficients are tabulated for 
different types of measurement schemas [7]. In addition, there is a simple formula 
discussed in the Appendix, that calculates averaging coefficients as function of the 
number and distribution of measurement levels, enabling the finesse analysis of the 
corrections. If no such detailed study is needed, one may employ other methods of 
accounting for disturbances, for example continuous least-square solution [8, 9, 2, 10], 
or shifted Legendre approximation [11, 12]. 



3. Distortion of registered trajectory due to finite speed of light 

3.1. Distortion as signal delay 

3.1.1. Delays in laser beam propagation Let for some moment r the vertical coordinate 
of the test body be z(t), and the distance from its origin to the beam splitter be b 
(fig. 2). Then the separation from the test body to the beam splitter is b + z{t), if the 




b 



Figure 2. Distortion of the registered trajectory resulting from signal 
delays, viewed cither as time delay At or coordinate 
progression Az. The case of the test body approaching 
the beam splitter. 

interferometer is positioned above, or b — z(t), if positioned below the test body. From 
here on, we combine both cases as b ± z(t), where the upper or lower sign corresponds 
to the beam splitter positioned above or below the test body, respectively. Due to finite 
speed of light, the information on the coordinate z(t) is delivered with some delay - 
in the moment t, for which 

f-r=*±iW, (11) 
c 
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where c is the speed of light. So, the coordinate z(r) is matched in data processing to 
the moment t, rather that r, causing the distortion of the trajectory. We proceed with 
establishing some relationships between t and r. 

3.1.2. Two scales of time We considered a laser beam as information carrier, traveled 
with some delay from the test body to the beam splitter. In this formulation, the body 
poses as the origin of the information, while the beam splitter poses as the information 
destination. The variables r and t can be considered as two scales of time. The motion 
of the test body is related to the scale r, whereas the registration of the motion is related 
to the scale t. 

The relationship between r and t is best presented as readings of moving clock 
compared to the readings of the clock resting with the observer, assuming that both 
clocks were synchronized before the motion started. Distortion of the scale r with 
respect to the scale t is caused by the fact that the moving clock readings are delivered 
to the observer with the variable delay (13). If the moving clock would also rest with 
the observer instead, the readings of both clocks coincided (dashed line in the fig. 2b). 

We use the following properties of the test body's trajectory in future analysis: 

(i) The [r , t] interval is much shorter than the trajectory's duration T: 

t - t < T. (12) 
This property enables replacement of r with t in (11): 

t-r= 6 ^, (13) 

c 

(ii) The distance traveled by the test body during the [r , t] interval is much shorter 
than the distance to the beam splitter: 

\z(t)-z(r)\ <6. (14) 

(iii) Change in the body's velocity during the [r , t] interval is minuscule compared to 
the change over the entire trajectory: 

i(r)-i(t)«0. (15) 

Simplifications coming from these properties are based on the fact that the timescales 
t and r differ in terms ~ 0(c~ r ) (11), so interchanging the scales in expressions having 
another c in denominator creates difference terms of only ~ 0(c~ 2 ), which are negligible 
as the test body's velocity is diminutive compared to the speed of light. 

3.1.3. The distortion as obtained from signal delay Two ways are possible to deduce 
the trajectory distortion from signal delays. 

Adjustment of time. We can proceed from the fact that the coordinate registered at 
the moment t corresponds to the true coordinate at the earlier moment r: 

z(t)=<r). (16) 
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Expressing r via t using (13), we get the following formula for the registered trajectory: 



Adjustment of coordinate. While the laser beam is traveling to the beam splitter, the 
test body keeps on moving along its path, so at the moment the beam reaches the beam 
splitter the body's true coordinate has progressed a bit compared to the registered one: 

z (t) = z(t) + Az. (18) 

As, according to (15) the body's velocity changes insignificantly over the path Az, we 
can use 

Az = z(t)(t - t). (19) 
Substituting this expression into (18) and applying (13), we find 

m = 4t)(lT i f)-b i( ^. (20) 

The above formula also follows from (17) by Taylor expansion , which makes both 
formulas equivalent for future analysis. The registered acceleration we find from (20) as 

g(t) = l{t) = z(t) =F Mt) c Z(t) - Z -^{b ± z(t)). (21) 
For uniformly accelerated motion 

z(t) = z + V t + g t 2 /2 (22) 
the registered acceleration, according to (21), is 

Q 

g(t)=goT—(V + g t). (23) 

c 

Remark. Equivalence of the formulas (17) and (20) allows the trajectory distortion to be 
obtained as either time delay or coordinate progression. Combining both approaches, 
however, will cause incorrect doubling of the distortion. 

3.2. Doppler distortion of the registered trajectory 

The reasoning we used to get the additional term of the acceleration in (23) was based 
on abstract considerations regarding signal delays. We did not rely on any specifics of 
a particular signal type, like passive or active, pulse or continuous, etc. Neither did we 
specify how the information is coded with the signal. Historically, however, analysis of 
the correction due to finite speed of light relied on the implemented physical phenomenon 
of interference of direct and reflected laser beams, manifesting itself as Doppler effect. 
The frequency of the light wave reflected from an object moving with the velocity V(r) 
after experiencing two Doppler shifts is given by [13] 

^^ c Mr <24) 
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where vq is the frequency before reflection. The frequency of the interference signal, 
neglecting terms ~ 0(V/c) 2 , will then be 

/ / \ s n V(t) 2u V(t) 2z/ 0ta/ , / V(t) 



c 



The velocity of the moving reflector can be found as 

V(r) = \v{t) = V(r) T (26) 

because Az/ = c, where A is the laser wavelength. For uniformly accelerated motion 
V(t) =V + g Twe get 

V(r) = V + go r^^±J^l. (27) 

c 

Acceleration found as dV/dr would be 

~ g (r)=g T^f(V + g r). (28) 

Additional acceleration term of this formula differs from that obtained earlier in (23). 
The reason is that we implicitly assumed the interference (25) to take place at the 
moment of reflection, while it's actually taking place with the delay defined by (13). To 
account for the delay, we need to convert the scale r into the scale t. Substituting (13) 
into (27), we get: 

As before, taking the derivative and dropping components ~ 0(V/c) 2 or less, we get 

g(t)=g oT ^(V + g t). (30) 

3.3. "Correction due to the Doppler effect" or "correction due to finite speed of light?" 

We are ready to highlight a misconception responsible for a significant part of 
discrepancies in the results, obtained by different authors. We first have to admit that 
both terms in the title of the current chapter are misnomers. Physical phenomena 
giving names to corrections are usually of a secondary nature, and their influence 
on the measurement result is relatively small. In this view, the Doppler effect is a 
primary phenomenon, making interferometric measurement of gravity at all possible. 
The fmiteness of the speed of light, in turn, is the reason for the Doppler effect to exist, 
so it's also a primary phenomenon. More accurate would be the correction terms like 
"non-linearity of the Doppler effect" and "non-uniform delays of the light signal". It's 
with these clarifications in mind, that we still use the old, historically settled terms in 
this paper. 

For a long time both terms were used as synonyms. In fact, the correction due to 
the Doppler effect accounts only for the | of the total correction due to finite speed of 
light, if, as shown above, no signal delay is taken into consideration. 
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If the Doppler effect is used for the correction, the "missing" | of the total value 
can be "recovered" by rendering the reflectors's velocity into the scale t. Indeed, if the 
actual velocity is 

V(t) = V + 9o t, (31) 
the observed velocity, because of (13), will be 

T/m T/,_ (* b±V t±g t 2 /2 \ 

V[t) = V + g lt I, (32) 

which is equivalent to the observed acceleration found as dV/dt: 

g(t) = goT^(V + g t), (33) 

representing exactly | of the additional component in (30). 

4. Compensating distortions from the finite speed of light 

4-1. Correcting the measured gravity 

By averaging the distorted acceleration (23) with the weighting function (7), we find, 
using (10), that the additional value measured by the gravimeter due to finite speed of 
light is 

Ag~ c = T — (V + goC 1 ). (34) 

c 

Taken with the opposite sign, the above expression yields the correction due to finite 
speed of light. We now find the correction for different measurement schemas. 

4-1.1. Two-level schema Substituting (B.2) into (34) and allowing the only sensible 
value of Vo = for this schema, we get 

Ag c = (35) 

c 

4-1.2. Three-level schema Substituting (B.3) into (34), we get 

Ag c = ±^(g (T 1 + T 2 ) + 3V ). (36) 

c 

4-1.3. Four-level schema Substituting (B.4) into (34), we get 

g f T 3 2 + T| — Tf + T 2 T 3 qt/ ^ 

&g c = ±— [go Tj, , rp 7f + 3 Vb . (37) 

c V I3 + 12- li J 

If T 3 -T 2 =T 1 , T 2 -T x = t, the correction will be 

Ag c = ± 3 ^ (g (27\ + r) + 2 V ) . (38) 
2c 
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EST schemas For levels equally spaced in time, the correction due to finite 
speed of light, according to (B.8), does not depend on the number of levels and equals 

A£ c = ±^(Vo + ^oT). (39) 
4-1.5. ESD schemas 



Accurate calculation of the correction. Because the averaging coefficient C\ in ESD 
case depends on the number of levels, the accurate correction value can be obtained 
by calculating the coefficient with the formula (B.13), and substituting it into (34), 
resulting in 

A 9c = ± 3 -^(v + 9 iG(T l ,T?)), (40) 



c V 6 

where the vector T iy % = 1...N is calculated using (B.14). 

Case N — > 00. For infinite number of levels, the calculations can be simplified using 
the C\ coefficient given by (B.15), which yields the following correction formula: 

A„ -+ 39o fv 1 ^(^o 3 T 3 + 45^T 2 K, + 108^T\/ 2 + 70 \/ 3 ) \ 
9c c V + 9 ° 7 (glT3 + l2glTiVo + 30g TV* + 20V 3 ) J " W 

Case N — > 00 and Vq — while calculating C\. Simplified formula (B.16) for C\ leads 
to the following correction: 

Ag c = ±^(v + ^g T), (42) 



7 

also given by Kuroda & Mio in [5]. It must be stressed here, that even though the above 
formula includes the term Vq, it does not completely account for the effects of initial 
velocity, because it was assumed to be zero when obtaining the averaging coefficient ^T. 
To see how accurate the corrections given by the simplified formulas (42) and (41) are, 
in the table 2 we find the corrections for some real instruments that use ESD schema. 
As follows from the table, the magnitude of the correction (42) by Kuroda & Mio [5] 
can for real instruments exceed the accurate value up to 1 /iGal. On the other hand, 
the approximation (41) for N — > 00 produces corrections accurate up to the 0.01 /iGal, 
if the number of levels is at least several hundreds. 

4-1.6. Symmetric rise-and-fall schema Because of the (B.17), the correction due to 
finite speed of light is negligible for the symmetric rise-and-fall schemas. 

4-2. Refining the trajectory model 

Due to the signal delays, the model of the uniformly accelerated motion 

^ = 20 + ^0^ + ^0^/2 + ^, (43) 
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Table 2. Corrections due to finite speed of light for some instruments using ESD 
schemas: 

(40) is the accurate correction. 

(41) is the approximation for N — > oo. 

(42) is the Kuroda & Mio's approximation [5]. 

schema parameters Ag c , fiGal, by formula: 
gravimeter setup . 







t/ m 


T, s 


N 


(42) 


(41) 


(40) 


JILA 


[14] 


0.30 


0.20 


45 


-13.94 


-13.35 


-13.11 


JILAg 


[15] 


0.20 


0.19 


170 


-12.41 


-11.95 


-11.87 


FG5 


[6] 


0.40 


0.17 


700 


-13.27 


-12.66 


-12.65 


A10 


[16] 


0.30 


0.09 


200 


- 7.89 


- 7.52 


- 7.51 


MPG-1 


[17] 


0.32 


0.20 


1635000 


-14.14 


-13.53 


-13.53 


MPG-2 


[181 


0.32 


0.16 


1100000 


-11.94 


-11.40 


-11.40 



where 

Tj are the measured time intervals, 

Si are corresponding measured distances, 

t{ are the measurement errors, 

is statistically inadequate, if the time intervals are related to the beam splitter. The 
least squares estimate (A. 6) of the parameter g in this case has a bias defined by (34). 
To bring the model adequate, the time intervals without delays should be used, i.e. the 
time intervals related to the reflector, not the beam splitter. This transition is inverse 
to what we did in the chapter 3.1.3, to find the correction (statistically: bias). The 
refined model is 

Si = z + V fi + g fi 2 /2 + e h (44) 

where 

f t = T l + (45) 
c 

The property (14) enables us to use the measured distances instead of the real ones. 
Time intervals corrected as above will result in unbiased least squares estimate of the 
parameter g , so no further correction is needed. The only caveat in this approach 
is to estimate the parameter b. One approach [6] would be to use z of (44) as b. 
Unfortunately, there is no way to arrange data acquisition and processing in absolute 
gravimeters so that parameter Zo would represent the distance from the beam splitter 
to the initial position of the reflector . Moreover, identifying b with z$ would render the 
model (44) non-linear with respect to zq. Another approach, suggested in [19], is to put 
6 = 0. The possibility to nullify b, as well as to assign it any arbitrary value follows 
from insensitivity of the estimate (A. 6) to any constant additive of time [20]: 

G(T t + A,S t ) = G(T t ,S t ), (46) 
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where A is some constant, which in our case equals b/c. The formula for the corrected 
time intervals would then be 

f i = T l ± ^. (47) 

The property (46) also holds true for the tree- and four-level models (A. 3) and (A. 4), 
so the simplification (47) can be used for these models as well, but not for the models 
(A.2) or (A.5). 



5. Corrections due to finite speed of light obtained by other authors 

5.1. Hammond and Faller [21] 

The correction reported by the authors for the three-level schema is 

A9c = "7 {1 90 (Tl + T2) + 2 v °) ■ (48) 

The discrepancy with the result (36) is partially explained by identifying the corrections 
due to the Doppler effect and due to finite speed of light, as discussed in the chapter 3.2. 
In that case the coefficient | should be |. 

5.2. Arnautov et al [22, 23, 24, 3] 

Using the phase progression, in the paper [22] the authors assessed the trajectory 
distortion similar to (36) and found the correction for the three-level schema same 
as (20). For the specific case of four-level schema we considered in chapter 4.1.3, the 
correction obtained by the authors [22, 23] is thrice less than in (38). 

In the paper [3], the authors applied the formula (A.5) to the four-level schema. 
The correction given by the authors corresponds to the coefficient C\—{Ti + T 2 + T 3 )/3. 
According to our estimates, the C\ in this case is defined by (B.5). 

5.3. Murata [2] 

The paper [2] deals with the EST schema, for which the correction we found is (39). 
The paper reports the following formula for the correction 

Ag c = \& T. (49) 

To get the correction, the Doppler effect (only) was applied to the two-level schema, 
making the result equal | of the correction for such a schema (35). 

5.4- Zumberge [25] 

The gravimeter described in [25] implements the ESD schema, but the applied correction 
agrees with that for the EST schema (39). To get the correction, the author averages 
acceleration like we do in (6), but with no weighting function. This averaging actually 
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corresponds to the uniform weighting function w(t) = ^, having the first averaging 
coefficient 

T 

i r t 

C 1 = T Jtdt = -, (50) 
o 

which coincidentally corresponds to the value for the EST schema (B.8), thus explaining 
the obtained result. % 

5.5. Kuroda and Mio [5] 

Considering both signal delay and Doppler effect, the authors found the interference 
phase, and integrated it to get the distorted trajectory. For the three-level schema 
(A. 3) the correction was found to be 

Ag c = -^(g (T 1 + T 2 ) + 2V ). (51) 

c 

For the similar schema described in the paper [26], the correction was determined as 
A 9c = ~j (9o (*i + *2 + h) + 3 V ) , (52) 

where t\ and ti are the start and the end of the T\ interval, tj, is the end of the T 2 
interval (fig. Al b). In the first formula, the coefficient of Vo should be 3. The second 
formula is valid when t\ coincides with the time origin. 

Asymptotic value of the correction the authors obtained for the EST schema is the 
same as in (39). For the ESD schema, the authors got 

Ag c = ±^(V + ^g T). (53) 
c 7 

The coefficient | of the above formula corresponds to zero initial velocity and infinite 
number of levels. For real shown in the chapter 4.1.5, the above correction 

may up to 1 /xGal exceed in absolute value the accurate one. 

5.6. Hanada [4] 

The paper evaluates the distortion of the registered trajectory based on signal delays. 
For a small increment of the observed time At, the corresponding distance Az is found. 
The total length and duration of the trajectory are determined by integration, and the 
acceleration is found using two-level schema (A. 2). One should expect this approach to 
yield the correction like (35). Indeed, the increment of the real time corresponding to 
the At increment of the observed time, can be found based on (13) as + 

At = (l + ^\ At. (54) 



This coincidence can also be explained by symmetry of the uniform function — sec formula B.9 
+ Following the author, we assume the trajectory of the test body to be z(t) — g r 2 /2. 
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The increment of the distance traveled in time At is 

Az = g TAr. (55) 

Expressing r in terms of t using (13) and substituting At from (54), we find the total 
distance registered in time T: 

*-/Va,-»/(. + £)( 1+ *)-.-!£ + ^. m 

o o 
The registered time is 

T T 

T — J At = J dt = T, (57) 
o o 

leading to the following registered acceleration 

g= 2 4 = 9o + 9 -^- (58) 
y c 

We've just very likely followed T. Tsubokawa to his unit value for the disagreement 
coefficient k in (1), as cited by Hanada on pages 99 and 134 of his paper. However, the 
value obtained by Hanada himself is |. The discrepancy is the result of the following 
assumptions. 

(i) The information-bearing signal is considered propagating from the beam splitter to 
the reflector, rather than vice versa, causing the opposite sign of the time distortion. 
The formula (3.6.22) of the paper (corresponding to our formula (54)) should have 
the opposite sign. 

(ii) Expressing the registered distance via the body's velocity, the time distortion is 
accounted for only in the time increment, but not in the time itself. The formula 
(3.6.24) of the paper treats the distance increment as g tAr, rather than g rAr, 
like in (55) 

(iii) The paper implements double accounting of signal delays, as mentioned at the end 
of the chapter 3.1.3. The paper's formula (3.6.23) treats the registered time as an 
integral of At, rather than that of At, like in (57) 

(iv) The factor 2 in the paper's formula (3.6.22) should be 1. 

While every assumption changes the result in its own way, their cumulative effect can 
be found with the formula 

T 

f(t + PD8$f) (l + D 9 -f)dt 

9 = ^"— T ^— , (59) 

}(l + PQD 9 -f) dt\ 

being just the two-level schema formula g = 2S/T 2 , in which the registered distance 
and time are expressed through the formulas (56) and (57) with the option to include 
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Table 3. The factors of the formula (60), emulating the influence of different 
assumptions made in the paper [4] on the correction due to finite speed of light. 



assumptions factors 



(i) 


(ii) 


(iii) 


(iv) 


D 
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R 




- 


- 
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1 
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l a 
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1 
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1 
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/ 
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-1 


— 




— 
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1 
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2 
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f 
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2 
3 
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1 
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1 
3 




■j 


</ 


■J 


1 
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1 
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4 
3 










-f 


1 
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-i 


✓ 
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-f 


1 
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-i 


✓ 




/ 




-f 
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1 


1 





/ 




✓ 


/ 


-f 


1 


1 


2 


1 


/ 


/ 






-i 








1 


2 
3 


/ 


/ 




/ 


-1 








2 


2 
3 


✓ 


/ 


✓ 




-1 





1 


1 


lb 

3 


/ 


✓ 


/ 


/ 


-1 





f 


2 


4 c 
3 



a The result of T. Tsubokawa. 

b The result of Hanada, no evaluation error. 

c The result of Hanada. 



or exclude the above assumptions. For example, the sign of time distortion (assumption 
(i)) is modeled by the factor D assuming values: 1 or -1 (see table 3). The above 
expression is equivalent to 

g = g {l + g -^L (^PD + ^D-RQD^y (60) 

So, the correction is influenced by the assumptions in the following way 

Ag c = -k 9 -^, (61) 

c 

where 

k=\pD + \D-RQD. (62) 

The table 3 shows possible values of D, P,Q, R , and the resulting values of k in formula 
(61), which cover the range from — | to | with the increment of |. 

5.7. Niebauer et al [6] 

The instrument described in the paper [6] implements refining of the trajectory 
approach, which, according to the chapter 4.2, gives accurate correction value. 



Correction due to finite speed of light in absolute gravimeters 



16 



Discussing the correction, the paper presents its approximate value as — §^-*, which 
corresponds to the Vq — case of the EST schema (39). As the instrument actually 
implements the ESD schema, more accurate approximation is given by the formula (42), 
i.e. — y^r- Two approximations differing 1.4 /iGal suggest some degree of uncertainty 
in the knowledge about the correction [27], which probably should have been reflected 
in the error budget, as it includes entries as low as 0.1 /xGal. 



5.8. Nagornyi [7] 

The approach suggested by Kuroda & Mio in [5] was followed with the weighting 
functions treatment, as in chapter 2, leading to the correction formula similar to (34). 
The results of Kuroda & Mio for EST and ESD schemas were confirmed as special cases 
of the formula (34). Though the paper [7] contained all necessary information about 
the coefficient C\ to write down the formula (41), no further analysis of the correction 
was done at the paper. 



5.9. Robertsson [11] 



Using the signal delays approach, the distortion component of the trajectory was found 
as in formula (20). Projecting the distortion component to the space spanned by the 
shifted Legendre basis, the correction for the case of a test mass falling away from the 
beam splitter was found as 

;i+%(a)/5) 



where 



where 



?72(A) 



3g V , 3Tg 2 



+ 



2c 



5A(A 2 - 3) 
7(3A 2 - 5) ' 



A = I 1 + 2V /g T 




for ESD schemas, 
for EST schemas. 



(63) 



(64) 



(65) 



If put together, the above equations yield the same correction formulas as (39) and (41). 



6. Conclusions 



In the paper we attempted to revise the correction due to finite speed of light in absolute 
gravimeters in order to understand and reconcile differences existing in the theory. Like 
other authors, we based our reasoning on two interrelated physical phenomena: delay 
in the electromagnetic wave propagation, and the Doppler effect. To achieve our goal, 
we had to make the following advancements in the implementation of the phenomena 
models: 

* We have fixed the obvious typo in the text. 
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(i) Introduction of two scales of time associated with moving and resting reflectors and 
deriving the set of rules allowing for easy transition between the scales; 

(ii) Establishing the equivalence between the observed trajectory distortions derived as 
either delays in registered time, or progressions of registered distance. 

The analysis we made leads to the following conclusions: 

(i) The Doppler effect in the scale of time associated with the test body is responsible 
for the | of the correction due to finite speed of light in absolute gravimeters. 
Transition of the test body's velocity to the time scale associated with the beam 
splitter, which is equivalent to introduction of optical signal delays, is necessary to 
get the full correction value.. 

(ii) For EST schemas implementing least-squares fitting of the three-parameter linear 
model (working formula (A.6)), the correction due to finite speed of light does not 
depend on the number of levels and always defined by the formula (39). 

(iii) For ESD schemas implementing the same model, the correction is a convoluted 
function of number of levels and initial velocity. Approximation by Kuroda & Mio 
(formula (42)) may be up to 1 yuGal off for real instruments. The approximation 
by formula (41) is 0.01 yuGal accurate, if the number of levels is at least several 
hundreds. 

(iv) If the finite speed of light is accounted for by correcting the measured time intervals, 
the parameter in time recalculations can be dropped for the following schemas: 

• three-level schema (A. 3), 

• four-level schema (A. 4), 

• any multi-level schema implementing least-squares fitting of the three- 
parameter linear model (A.6). 

(v) Discrepancies in the corrections obtained by different authors are caused by several 
reasons, including confusion between corrections due to the Doppler effect and due 
to finite speed of light, use of simplified trajectory models, and math errors. 

Appendix A. Weighting functions of absolute gravimeters 

Appendix A.l. Working formulas of absolute gravimeters 

Weighting function of absolute gravimeters are based on their working formulas that 
map measured from the common start intervals of time {Ti\f =l and distance {Si}f =1 to 
the measured gravity acceleration}]: 

g = G(Ti,Si). (A.l) 

The term "working formula" we use here is narrower than the term "measurement 
equation" , which usually includes all the calculations leading to the measurement result 

U Because intervals are measured between levels, the number of levels is always one more that the 
number of intervals. 
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[27]. It's essential that working formulas are linear in measured distances. In the paper, 
the following working formulas are considered: 

• Two-level schema (fig. Ala, N = 1) 

g = G(T,S) = 2 ^. 

• Three-level schema (fig. Alb, N = 2) 

5*2 _ 5*1 

T> ' T. ) 7. / 

3) 



g = G(T t , S t ) = 
• Four-level schema (fig. Ale, iV 
g = G(T t , S % ) = 



(A.2) 



(A.3) 



S?j — So 



Si 
Ti 



T 3 + T 2 - T\ 



T 3 -T 2 

• Multi-level schema with least-squares fitting of two-parameter model Si 
9ot 2 /2 

1 En E^Si 

ETf ET?S t 



g = G(Ti, Si) = 2 



En E T i 

et! En 



Multi-level schema with least-squares fitting of three-parameter model Si 
V t + g t 2 /2 



(A.4) 
V Q t + 

(A.5) 
= z Q + 



g = G(T l , Si) = 2 



N ETi ESi 

ETi En ETiSi 

En En En Si 



n ETi En 
E T i En En 
En En En 



.(A.6) 



Appendix A.2. Determination of weighting functions 

To find the weighting function of a gravimeter with working formula G (Tj , Si) linear in 

distance intervals Si, one should replace Si with piecewise linear functions hi(t) defined 
as [7] (fig. A2) : 

f Ti-t 0<t<Ti, 

hi(t) = < ~ ~ A.7 

that means 

w(t) — G(Ti, hi(t)) . (A.8) 

Weighting function for the two-level schema (A.2), will therefore be generated by a 
single component like (A.7): 

w(t) = 2 7 ^. (A.9) 

This weighting function is shown on the fig. Ala. For the three-level schema (A.2), 
the weighting function is generated, according to its working formula, by a linear 
combination of two scaled components like (A.7). This process is shown in detail on the 
fig. A3 and results in the triangle-shaped weighting function (fig. Alb): 
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T 2 T 3 



Figure Al. Measurement schemas of absolute gravimeters and corresponding weighting functions: 
a: two-level schema; 
b: three-level schema; 

c: four-level schema, the weighting function depends on the working formula: 

1 - working formula (A. 5), 

2 - working formula (A. 6), 

3 - working formula (A. 4). 



Ti 




Figure A2. The function hi(t) (A. 7) that generates gravimeter's weighting function, 
when substituted for the Si interval to the gravimeter's working formula 
G (Ti , Si). 



2t 
T\T 2 



w{t) = < 



2(T 2 -t) 
I T 2 (T 2 - 



0<t<T u 



T x < t < T 2 



(A.10) 
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Ti 
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V T 2 Ti ^ 



T 2 -T! 




Figure A3. Formation of the weighting function for the three-level schema 
For the four-level schema (A. 4) the weighing function is trapeze-shaped (fig. Alc-3) 
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w(t) 



T^Ts + Tz-n) 
2 

T 3 + T 2 - Ti 

2(T 3 -t) 



0<t<T u 
Ti < t < T 2 , 
T 2 <t< T 3 . 



(A.11) 



i (n-TzXn+Ti-T!) 

In the same way, weighting functions of multi-level schemas are linear combinations of 
N components like (A. 7), resulting in piecewise linear weighting functions with N — 1 
segments. As shown on the fig. A4, the j-th segment of the weighting function is created 



L 3-l 



T 2 




Figure A4. For the formula (A. 13): on the time interval Tj. . . Tj + \ only those hi(t) 
components are non-zero, for which i > j. 



only by those hi(t) components, for which i > j. Than is, 
Wj(t) = G(Ti, Sji), 



where 



Si 



XI 





Ti — t 



1,3, 



j + 1 , N. 



(A.12) 



(A.13) 



Weighting functions obtained with the formulas (A.12) and (A.13) are shown on the 
fig. A5. As the number of levels increases, the functions converge to a limiting shape. 
For the EST schema with working formula (A. 6), the limiting shape is defined by 



t 2 t 3 t A 
w (t) = 30— -60— + 30 — . 

V / rp3 rp^ rp^ 



(A.14) 
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For the ESD schemas with working formula (A. 6), the limiting shape depends on the 
initial velocity and is defined by the equation: 

60£ 4 120t 3 60t 2 \ /A x 

-jr, fT- + -rr) (A.15) 

gltT 2 + 3g 2 T (2t + T) V + 6 g (t + 2T) V 2 + 10 V$ 



gl T 3 + 12 g 2 T 2 V + 30 g T V 2 + 20 V 3 

Weighting functions for the ESD schema approach those for the EST schema as Vo — > oo, 
but even for small initial velocities they are significantly deviating from the case Vq = 
(fig. A6). A method of finding limiting shapes can be found in [7]. 

The working formulas (A. 5) and (A. 6) can also be applied to the schemas with 
three and four levels. For three levels, the weighting function is still the same triangle 
(fig. Alb), as there is only one parabola passing through three points. For four levels, 
the weighting function significantly depends on the working formula (fig. Ale). 
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Appendix B. Averaging coefficients of absolute gravimeters 

Based on the generic formula for weighting functions (A. 8), the averaging coefficients 
can be found as [7] 

Cn ~ (n + l)(n + 2)- 

For the correction due to finite speed of light, only the C\ coefficient is needed. According 
to (B.l), this coefficient is 

• For the two- level schema (A. 2): 

Ci = |. (B.2) 

• For the three-level schema (A. 3): 

Ci = (B.3) 

• For the four level schema (A. 4): 

Tj + Ti-T? + T 2 T 3 

Cl = m+T 2 - Tl) ■ (B - 4) 
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• For the schema with four levels and working formula (A. 5): 

= 1 T 1 2 t|(T 1 -T 2 ) 2 (T 1 +T 2 )+T 1 2 t2(T 1 -T3) 2 (T 1 +T3)+t|t|(T 2 -T3) 2 (T 2 +T 3 ) /g g\ 

1 3 T 2 T 2 (7\ -T 2 ) 2 +T 2 T 2 (T! -T3) 2 +T 2 T 2 (T 2 -T 3 ) 2 V ' / 

• EST schemas. If levels equally spaced in time, the i-th time interval is 

Tj = T ^— — (B.6) 

where T is the duration of the measurement interval. Substitution of (B.6) into 
(B.l) creates to the following sums: 

i=l v 7 

For different ra'a, the values of the sums are available in references [28]. Substituted 
to the formulas for C 1 (B.l, A. 6), the values result in 

Ci = \ (B.8) 

Another way to prove this fact would be to observe, that regardless of the number 
of the levels, weighting functions for the EST schemas are always symmetric with 
regard to the middle point (fig. A5), so 

T 



I 



t w(t) dt = J (t-l^) w ( t ) dt + \ J w ( t ) dt = \- ( B - 9 ) 





The first integral equals to because of the symmetry, the second integral equals 
1 sue to the unit property (5). 

• ESD schemas. If iV levels are equally spaced along the total path S, the i-th level 
can be found as 

s > = s jfri- ( R1 °) 

The time interval, corresponding to the Si can be found using 

S i = V T i + g T?/2. (B.ll) 

This equations doesn't include initial displacement, because for calculations the 
origin of the time axis can be co-located with the first measured levelff. Using the 
length of the trajectory 

S = V T + g T 2 /2, (B.12) 

where T = T N is the total measurement time, we get 

-V + A 2 + 2g (V T + g Ty2) g 

T,= * . (B.13) 

9o 

ff This shift of the time origin, however, does not equate the models Si = V^Ti + g Tf /2 and 
Si = zq + VoTi + goTf / '2 when used for the least-squares fitting of the trajectory. These models 
have different weighting functions in terms of (4) and so they average the disturbances in the test body 
motion differently [7]. 
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This Ti vector we need to substitute into the equation (B.l), which for n — 1 
becomes 

C^G^T*), (B.14) 

to obtain the C\ coefficient, which turns out to be a complicated function of g , 
Vo, N, and T. If AT — ^ oo, the C\ can be found using the limit weighting function. 
Substituting (A. 15) into (10), we get 

T (4 gl T 3 + 45 gj T 2 V + 108 g T Vj + 70 Vg) 
° 1_ 7(^T3 + 12^T2\/ + 30^T\/ 2 + 20 \/ 3 ) ' 1 ' °> 

For the case V = 0, the above expression simplifies to 

Ci = ^T. (B.16) 

In the chapter 4.1.5 we investigate how accurate the two above formula are when 
applied to real instruments. 

For the symmetric rise-and-fall measurement schemas 

d = 0, (B.17) 
which follows from the evenness of weighting functions for such schemas [7]. 
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